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Abstract: In 'super singular' scattering the potential g 2 UA(f) involves a variable 
nonlinear parameter A upon the increase of which the potential also increases 
beyond all limits everywhere off the origin and develops a uniquely high level of 
singularity in the origin. The problem of singular scattering is shown here to be 
solvable by iteration in terms of a smooth version of the semiclassical approach to 
quantum mechanics. Smoothness is achieved by working with a pair of centrifugal 
strengths within each channel. In both of the exponential and trigonometric re- 
gions, integral equations are set up the solutions of which when matched smoothly 
may recover the exact scattering wave function. The conditions for convergence 
of the iterations involved are derived for both fixed and increasing parameters. 
In getting regular scattering solutions, the proposed procedure is, in fact, supple- 
mentary to the Born series by widening its scope and extending applicability from 
nonsingular to singular potentials and from fixed to asymptotically increasing, 
linear and nonlinear, dynamical parameters. 
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I. Introduction 



We are going to consider scattering by singular potentials g 2 UA(r) at fixed 
as well as asymptotical values of the linear and nonlinear parameters g 2 and A, 
respectively. The pioneering work in treating singular scattering by exact means 
has been done by Calogero 1 in terms of the phase approach to quantum particle 
dynamics. In particular, he calculated the high energy limit of scattering by inverse 
power potentials 2 . A further correct approach is due to Froemann and Thylwe 3 , 
who analytically proved the exactness of the first order WKB (Wentzel-Kramers- 
Brillouin) approximation for the case of pure inverse power potentials in the short 
wavelength limit. Esposito 4 worked in singular scattering by the wave function 
polydromy method to get a recursion formula between solutions with potentials 
of different inverse powers in the radial distance. Notice all the above approaches 
relied on nonlinear first order or linear second order differential equations solvable 
analytically, in terms of special functions, for certain simple potential shapes. 

Nevertheless, it is the method of integral equations that is the first rank 
candidate for solving two body scattering problems involving a general shape and 
stage of the potential singularity. The pertinent Volterra type equations can be, 
perhaps, best classified in terms of the reference scattering problems implied , as 
follows. 

The most evident reference basis seems to be the unperturbed case g 2 = 
taken at the physical energy k 2 and the physical orbital angular momentum /. The 
relevant integral equations set up for the regular wave function u + (r) work in the 
case of nonsingular potentials, exclusively (Born series). It is true that for finding 
irregular solutions u~ (r) , there are analogous integral equations available, for cases 
of nonsingular and singular potentials alike, as shown by Newton 5 . However, all the 
above mentioned procedures can only be employed at fixed sets of the scattering 
parameters. Namely, the residual potential contained is there invariably given by 
the expression A(r) = g 2 UA(r) , whence the iterated series blows up term-by- 

2 



term in the limit g 2 — > oo (strong coupling). The same is true of the supersingular 
limit. Namely, the form factor Ua{t), besides being singular for r — > at fixed 
parameters, increases, by definition, beyond all limits at each fixed point off the 
origin for A —> oo. One has therefore to seek alternative approaches outside of the 
Born series for solving these asymptotical cases. 

An interesting choice of reference basis, as proposed by Newton 5 , offers itself 
in the exceptional cases of the solution being known in closed form at zero energy 
and zero orbital angular momentum. The residual potential is then k 2 — 1(1 + 
l)/r 2 . It is worthwhile to note that the iteration of the integral equation converges 
exclusively for singular potentials. Moreover, convergence obviously holds in both 
limits g 2 — > oo and A — > oo. 

Quite wide is the scope of applicability of the third type of reference prob- 
lems which we are going to consider in detail. It is, in fact, the semiclassical ap- 
proximation. The Langer version of this method, see e.g. Newton's monograph 5 , is 
specified by the reference centrifugal strength (/ + |) 2 . Owing to this very choice, 
the model wave function reproduces the exact regular solution of quantum me- 
chanics (QM) near the origin in the case of nonsingular potentials. As to singular 
potentials, the small distance QM wave function and its WKB approximation do 
coincide, even independently of the selection of the reference centrifugal strength. 
The problem of a possible convergent iteration and that of asymptotical values 
taken by the potential parameters had remained thereby still open. Recently, 
Dolinszky 6 proposed a smooth version of the semiclassical approach which is free 
of the highly inconvenient turning point singularity inherent in both the Langer's 
and the standard WKB approximation. The new procedure has been applied for 
developing convergent series expansion of the wave function describing scattering 
by singular potentials at fixed 7 as well as increasing 8 linear dynamical parameters. 
In the present paper an improved and generalized version of that approach will 
be proposed, which also includes the asymptotical case of increasing nonlinear 
parameters present in the Schroedinger equation. 
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II. Physical and reference problem 

The physical problem to be solved first is two-body scattering by a central 
singular potential say, g 2 U A (r) at the energy k 2 in the channel of index /, with 
A representing a nonlinear parameter. The potential should fulfill the following 
requirements: 

[r 4 U A (r)] -> oo, [r -> 0]; U A (r) > 0, U' A {r) < 0, [r > 0]; 

(2.1) 

[r 3 U A (r)] -> 0, [r -> oo]; [/^(r) -> oo, [r > 0, A -> oo]. 

In order to construct a suitable reference system for treating the problem, we 
introduce a triad (A 2 , A 2 , A 2 ) of auxiliary centrifugal strengths, subject to the re- 
strictions 

A 2 (/) = ^[A 2 (Z)+A 2 (Z)], A 2 (/)>A 2 (/). (2.2) 

The constants A 2 and A 2 are for the present freely chosen. The concept 'matching 
distance' is defined as the positive root of the 'master equation' 

k 2 R 2 - g 2 R 2 U A (R) - A 2 = 0; R = R(k 2 , g 2 , A 2 ; A) > 0. (2.3) 

The radial distance r will be in general substituted for by the dimensionless radial 
coordinate t. We put 

[r>0]. (2.4) 
Different regions of the space are distinguished as follows: 
t < 1, region e, (exponential region); 

£ > 1, region r, (trigonometric region); (2-5) 
t = 1, (matching point). 

The scattering process is governed by the radial Schroedinger equation 

+ fc 2 tf 2 - g 2 R 2 U A (Rt) - l ^^-}u ± (t) = 0. (2.6) 



A regular- irregular pair u^ 1 (t) of its solutions behaves near the origin as 

w ± (t) -> [g 2 R 2 U A (Rt)}~* exp[±#.R J dt'U A (Rt')*], [*->0]. (2.7) 

The reference problem will be a special smooth version of the zero-order 
semi-classical approximation. The entire argument hinges upon a pair of wavenum- 
ber function squares for the exponential and trigonometric regions, respectively, 
as follows: 

K 2 (t) = K 2 (t) = -k 2 + g 2 U A (Rt) + [t < 1]; 

5 ( 2 - 8 ) 

= K 2 (t) = k 2 -g 2 U A (Rt)-^, [t>l]. 

These definitions imply a number of properties concerning behavior around the 
matching point. In particular, one extracts from the master equation (2.3) that 

[Kl{t)] t=1 = [K 2 (t)] t=v [K 2 (t)] t=1 = *L^. (2.9) 

Overall properties also follow from Eq.(2.8). Indeed, 

K 2 (t)>0, mm{K 2 (t)}= [K 2 (t)] t _ v [0 < t); 

(2.10) 

K 2 (ti) > K 2 (t 2 ) if £i<£ 2 ; K 2 T {U) < K 2 T {U) if t 3 <t 4 . 

These statements issue from the restrictions (2.1) imposed on U A (r). Small and 
large distance behavior can also be thence extracted such as 

K 2 (t) ^oo if t^O; K 2 (t)^k 2 if t -> oo. (2.11) 

For the derivatives of i^ 2 (t) we shall apply the notation: 

4 S) (t) = ^ 2(t) d ; U , [ a = l,2; 7 = 6,r]. (2.12) 
Hence one obtains by Eqs.(2.1)-(2.2) and (2.8) 

[£>?>(*) - DP(t)] f=1 = 2 4 {2A 2 - g 2 R 3 U' A (R)} > 0. (2.13) 

5 



A regular- irregular pair of reference wave functions w^(t) is denned as 

wf(t) = nt(t)exp[±ut(l,t)], [t < 1]; (2.14) 

w$(t) = ?? r (t)[C ± cosw T (l,t)) + S ± sinw T (l,f)], [t > 1], (2.15) 
where we introduced the 'amplitude function' and the 'phase function' such as 

Vj(t)= (^y) 4 ^ [7 = e,r]; (2.16) 

tv y {t 1 ,t 2 ) = R [ 2 dt'lK^t')], [7 = e,r]. (2.17) 

As to the constants C ± and ^ , these parameters can be, for the moment, 
freely chosen. Out of them, C + and S + will be specified later by smoothness 
requirements. The choice C~ and S~ will, in turn, remain once for all free but 
the only restriction 

C + S-S + C-^0, (2.18) 

which warrants independence of the functions w^(t) and w~(t)). Yet, there exists 
a sophisticated definition of the constants C~ and S~ in terms of C + and S + , 
namely the one implied in the identity 

r* dt' 

^(ijE^iKl-/^}. [7 = 6,r]. (2.19) 

This relationship is automatically satisfied in the region e by the definition (2.14). 
As regards the point t = 1 ± 0, Eq.(2.19) guarantees there smooth matching of 
w~(t) to w~(t) whenever w+(t) matches there w+(t) smoothly. The particular 
choice (2.19) is quite irrelevant from the viewpoint of the present argument. Nev- 
ertheless, it justifies the use of the same superscripts (~) over w T ~(t) and w e ~(t). 

The reference wave functions w^(t), [7 = e, r] , of Eqs.(2.14)-(2.15) solve 
the pair (e, r) of differential equations 

+ k 2 R 2 - W±(t) - l -^-}w±(t) = 0, [ 7 = e, r]. (2.20) 
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Notice the differential equation (2.20 is common for w+(t) and w~(t), in the 
exponential and trigonometric regions alike. Namely, calculation yields for the 
reference potential 

W±(t) = W^t), [ 7 = 6,r]. (2.21) 

with the notation 

W 7 (t) = g 2 R 2 U A (Rt) + A 7 (t), [ 7 = e,r]. (2.22) 

The expressions of the residual potential A 7 (£) introduced here are extracted from 
Eqs.(2.14)-(2.15) and (2.20) as 

A e (t) = —[D^Ht)] 2 + \Df\t) - X " {t) " * ( * + 1} , [f<l], (2.23) 

K(t) = -^[^(t)] 2 + ^ 2) (t) - A '~^ + 1) , [f > 1]. (2.24) 

Recall that K 2 {t) of Eq.(2.8) is by Eqs. (2.9) continuous at t = 1 but, due 
to (2.12)-(2.13), not smooth. Therefore, the residual potential develops at the 
matching point a jump. Indeed, 

A T (t + 0)-A e (t-0)^0. (2.25) 

For simplicity, see expression (2.24), we shall work hence forward with the cen- 
trifugal strengths 

\l = (l + \)\ \ 2 T = 1(1 + 1). (2.26) 

Observe that this choice is compatible with the inequality contained in the postu- 
lates (2.2). 

III. A pair of convergent expansions 

A comparison of the Schroedinger equations (2.6), set up for the exact wave 
function w + (r), and Eq.(2.20), solved by the semiclassical wave function w + (r), 
suggests construction of a pair of integral equations, 

vt(t)=wt(t)+ /Va £ (OG+(M>+(0, [*<1], (3-1) 
Jo 
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vtit) = «,+ (*) + J*dt>A T (t')Gt(t,t')v+(t>), [t > 1]. (3.2) 

The solutions v+(t) and v^(t), if exist, are solutions of the exact Schroedinger 
equation within the respective regions. In particular, vf(t) is uniquely defined 
by Eq.(3.1) and furnishes a regular solution of the differential equation (2.6) in 
the exponential region. The solution v+(t) of Eq.(3.2), in turn, while solving the 
Schroedinger equation in the trigonometric region, still involves two free constants, 
C + and S + , as noticed following Eq.(2.15). Smoothness requirement for the overall 
regular solution v + (t) at t = 1 is just sufficient to unequivocally specify these 
coefficients. As to the notation in (3.1)-(3.2), the residual potentials A 7 (t) have 
been defined by (2.22)-(2.24). The resolvents involved in the integral equations 
are formally given as 

G+(t,t') = -L[ w +(t)w-(t') - w-(t)w+(t% [ 7 = e,r], (3.3) 

where the Wronskians contained are by Eq.(2.21) constant and read in general 

ci+ = W^{w+(t);w~(t)} = const., [ 7 = e,r]. (3.4a) 
In particular, one obtains after some calculations 

d+ = -2kR, d+ = kR(C+S- - C~S+). (3.46) 
The general expressions (3.3) can be recast in terms of the local wave numbers as 
Gt(t,t') = 2 sinh ^(M0] [o <*><*<!], (3.5) 

R[K*{t)K*(V)]V ~ 

Irrelevance of any special choice of the basis wf{t) for inclusion in the formula (3.3) 
manifests itself in the absence of the constants C ± , 5 ,± from the last expression. 
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The solution of the integral equations (3.1)-(3.2) rests upon the recursion 
schemes 

wt n {t)= ftfAcWG+Ww+^tf), [n= 1,2,3..]; «,+ (*) = «,+ (*), (3.7) 
Jo 

u;+ n (t) = |di'A T (iOG+(t,fVj m -i(f) I [m=l,2..]; w^(t) = w+{t). (3.8) 

A necessary condition for getting the solution sought for by iteration is convergence 
of each of the infinite series 



'(*) = » [t<l]\ (3.9a) 



n=0 



v+(t) = w tm(t), [*>!]• (3-96) 

m=0 

The overall solution v + (t) = {v+(t); v+(t)} should be throughout smooth. This 
requirement is realized off the matching point spontaneously. At the matching 
point, it is the zero order term of the series (3.9b) that exclusively contributes to 
both the solution v+(t) and its first derivative. If the series (3.9a) is in the region 
e convergent, the smoothness postulate for t = 1 can be recast in a simple form as 

[wt(t)]t=i = k + (f)]t=i; [wt(t)']t=i = k + (f)W. (3-10) 

These conditions simultaneously fix the trigonometric constants C + and S + thus 
completing the solution in the region r. The convergence proof for the above series 
becomes more transparent by using the following notation in both regions 7 = e, r 

*® = * w = 1 - (3 ' n) 



(t,,ta)= r it'\ Pl (t')\. (3.i3) 
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. In the exponential region, the formula (3.7) can thus be rewritten as 

q+(t) = /V^OIl-exp^MmVi^')- (3-14) 
Jo 

Hence we get by Eq.(2.17) the inequality 

\qtn(t)\ < f <tf)ptf [t < 1]. (3.15) 

JO 

Iteration yields then by Eq.(3.13) 

\Qen(t)\ < l[P e (0,f)] n , [t<l], (3.16) 
whence one extracts by analysis 

oo 

l«4(*)l < \w £ (t)\exp[P e (0,t)], [t < 1]. (3.17) 

The series v*(t) of Eq. (3.9a) is thus absolutely convergent if and only if the integral 
P e (0,t) exists and is bounded in t = (0,1). As to the trigonometric region, the 
resolvent formula (3.6) combines with the notation (3.12) to an equivalent form of 
the recursion relationship (3.8). Accordingly, we get 

t i 

w+n(t)= f d^(Osinu; T (M0^^™+ m -i(O, [*>!]• (3-18) 

The monotonicity relationship (2.10) implies then the inequality 

\w+ m (t)\<^ dt'lp^Ww^itX [t>l,m= 1,2,3...]. (3.19) 

Recall now the definition (2.15) and conclude by the zero order, m = 0, identity 
in the relationships (3.8) that 



\w 



Ut)\<(kR)H\C + \ + \S + \), [t>l\. (3.20) 
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Iteration furnishes then by the inequalities (3.19)-(3.20) 

\w+ m (t)\ < (kR)H\C+\ + \S+\)^-[P T (l,t)] m , [t>l;m= 1,2,3...]. (3.21) 

Summation over m yields thus 

00 

\< m (t)\ < (kR)H\C + \ + \S + \)e W P T (l,t), [t>l]. (3.22) 

m=0 

So the series v^(t) of the definition (3.9b) is absolutely convergent whenever the 
integral P T (l,t) of Eq.(3.13) does exist. 

Suppose the regional existence conditions 

P e (0,t)<oo, [£<1], P T (l,t)<oo, [t>l], (3.23) 

and, in addition, the smoothness postulate (3.10) are satisfied. If so, then solutions 
of the integral equations (3.1)-(3.2) together recover the regular solution of the 
differential equation (2.6). In particular, 

u+(t)=v+(t), [t<l]; u + (t)=v+(t), [t>l]. (3.24) 

Virtually, one always works , instead of the infinite expansions (3.9a)-(3.9b), with 
cut-off series such as 

N 

v^ N \t) = J2<n(t), [t<l]; (3.25) 

n=0 
M 

v^ NM) (t)=J2<i N) (t), [t>l\. (3.26) 

m=0 

Smooth matching of these functions at t = 1 in terms of the relevant analog of 
Eqs. (3.10) fixes the trigonometric coefficients C + and S + in terms of the cut- 
off 'length' iV in the exponential region but independently of its pair M in the 
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trigonometric region. These constants therefore should carry the superscripts (N) 
only. The postulates (3.10), rewritten for the cut-off approach, read 

2*(kR) 1 *C + W = [vt {N) (t)] t=u (3.27) 

-2-i[Di 1 \t)] t=1 (kR)^C+^ = [v+W(ty] t=1 . (3.28) 

Notice the cut-off approach must not be used unless the conditions of convergence 
are fulfilled in both regions e and r. 

IV. The supersingularity limit 

The term 'supersingularity' implies in our terminology a singular potential 
which is subject to the requirements (2.1). This involves, among others, that 
it contains a nonlinear parameter A which may increase beyond all limits. The 
relevant asymptotical form of the master equation (2.3) reads 

k 2 - ^U A (R A ) - |J- - 0, [A^oo]. (4.1) 

Bounded values such as Ra — > and Ra — > const, are by the properties (2.1) 
obviously excluded from the large-A solutions of the Eq.(4.1). One is thus left 
with the only possibility that Ra — > oo for A — > oo. The master equation itself 
reduces therefore in the supersingularity limit to 

U A (R A )^^, [A^oo]. (4.2) 

It is perhaps worth recalling that U a (t) ^ if A = fixed, r — > oo while Ua{t) — > 
oo if r = fixed, A — > oo, as implied in the set of postulates (2.1). In Eq. 
(4.2), in turn, the singularity parameter A and the matching radius Ra vary 
simultaneously and both increase to infinity so as to keep the left hand side of the 
equation constant. 

Four classes of strongly singular potentials will be introduced with each 
potential being specified by a variable core parameter A and a fixed tail parameter 
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B. The dimensionless form factor is in each case a product of an, for r — > 0, 
exponentially or powerlaw increasing core factor and an, for r — > oo, exponentially 
or powerlike decreasing tail factor. Concerning the subscripts, we use an obvious 
notation when writing 

tt ( \ ( a P l l 3r \ 
Uaf3{r) = exp( ), 

r P2 
t/ a/3 (r) = (l + ^) a exp(^), 

^(r) = exp(^)(^, ^ 

V*{r) = {l + r +n-^-)\ 

r r2 + r 

[0 < r, < a,/3,r 1 ,r2,pi,p2] a > 4; b > 3]. 

In the last section, we established general criteria for the convergence of the 
expansions (3.9a)- (3.9b). In the present section, fulfillment of those requirements 
will be checked for increasing nonlinear parameters. 

In the supersingularity limit A — > oo,the asymptotical form of the master 
equation (4.2) is explicitly solvable in cases of exponential tail for both types of 
core singularity. In fact, one concludes from the definitions (4.3) that 

R a p -> (piP2)H^) 1/2 , [a->oo]; 

R a p -> (riP2)Hg)K [a->oo]. 

In the powerlaw tail cases, in turn, only implicit expressions can be extracted from 
the formulas (4. 2)- (4. 3). Indeed, one gets 



R a b -> r 2exp(— — ), [a-^ooj, 
ari 

«ab -»• r 2 exp(— — ), [a^ooj. 
on ab 



(4.5) 



The matching distance increases in all of our examples slower than the respective 
variable parameter A = a or a. Indeed, one finds by the relationships (4.4)- (4.5) 
that 

^^^0, [A^ oo; A = a,a;B = p,b}. (4.6) 
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The notation A — > oo and Rab oo compete in representing the super- 
singular limit by a single symbol. We prefer the use of the latter alternative 
where possible. Accordingly, we shall throughout eliminate from the formulas the 
parameter A in terms of Rab- 

It is also obvious by analysis that the present scattering formalism should, 
in the supersingularity limit, become, mutatis mutandis, identical for exponen- 
tially and power law increasing cores. The equivalence is realized at the following 
correspondence of potential parameters: 

r\a — > p\a, [a, a — > oo, B = f3, b]. (4-7) 

In the limit considered, one has thus to treat, out of the four cases in (4.3), only 
two essentially different ones, namely (A,/3) and (A, b). We now rewrite the wave 
number squares of Eqs.(2.8), in terms of the variable t of the definition (2.4), 
for Rab oo. In the cases (A, f3), this transformation is, owing to the explicit 
relationships (4.4), straightforward. Indeed, one obtains by the definitions (4.3) 
in the exponential region 

A'^W - ^^Ot - «)1 + - 1} ' (4.8) 

[t < 1, Ra/3 — > oo, A = a, a]. 

A little bit more complicated is the incorporation of formula (4.5) into Eq.(2.8) in 
power tail cases for which one gets still in the region e 

^w-^(^]-^-i + F |^}, (49) 

[t < 1, RAb — > oo, A = a, a]. 

In the trigonometric region, distinction should be made between S-wave and higher 
partial waves. Indeed, one extracts from the definitions (2.8) and (4.3) for the 
potential classes A = a and a equally, that 



K 2 TA( 5(t) - k 2 {l - -f^exp[-^(t - i)]}, [t > 1,1 = 0,R AP - oo], (4.10) 

ft Tq P2 t 
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K 2 TAb {t) - k 2 {l - [^(J^W ±}, [t > 1, 1 = 0, iU* - oo], (4.11) 

K 2 AB (t)^k 2 {l- fc2 ^| t2 >. [t>l,i>0,i2AB-oo]. (4.12) 

Notice for the highest values of the parameter A the wave number function becomes 
independent of the potential: 

K 2 AB (t)^k 2 , [*>l,Z>0,iUs->oo]. (4.13) 

The quantities Dj S \t) of the definition (2.12 will be calculated below from 
the set of supersingularity expressions (4.8)-(4.12). In the exponential region one 
finds for both cases A = a, a in the limit Rab — > oo 

B&(«)--^(p + D. B&W-^^ + D", [«>1,J>0], (4.14) 

O e %W--^ln(^), D'iW - ^H^f )] 2 , [<>1,(>0]. (445) 

As regards the trigonometric region, the formulae are again sensitive to the orbital 
angular momentum. Indeed, for the potential classes A = a, a alike, one gets 

^>)-^f(p + l), ^-(^ + 1)', [«>1,I = 0], (4.16) 

D TM>(f) -> JilH—)?, - -SH— )] 2 ' lt> 1,1 = 0], (4.17) 

^!1bW-t4^73' ^ 2 LW-"#t, [t>l,l>0]. (4.18) 

The residual potentials (2.23)-(2.24) are quadratic and linear expressions of 
the quantities (4.14)-(4.18). They read in the exponential region, for exponential 
and powerlaw cores in like manner, 

Ae^(f)->~(^) 2 (4 + l) 2 , [f<l,i^->oo], (4.19) 

ID P2 t 

15 



A e ^(t)^-^[ln(^)] 2 ^, [t<l,R Ab ^oc]. (4.20) 

The analogous formulae of the trigonometric region are sensitive to the orbital 
angular momentum as 

A TA/3 (t)--^(^) 2 (^ + l) 2 , [t> 1,1 = O,Ra0^oo], (4.21) 

ID P2 £ 

A r ^(t)--^[ln(^)] 2 ^, [t > 1,1 = (4.22) 
3A 2 

&rAB(t) - 2k 2 R 2 t ^ [* > 1, i > 0, Rab - oo, B = p, b]. (4.23) 

The next step towards checking expansions (3.9) for the realization of con- 
vergence criteria is calculation of the quantities p-y{t) introduced by Eq.(3.12). 
The asymptotical expressions (4.8)-(4.13) combine with the ones of (4.19)-(4.23) 
to yield, first for the exponential region, 

P^^-^^^^ + l^eM-^^-t)], [t<l,R A ^oc], (4.24) 
lb gp 2 t z p 2 t 



^(t)---[ln(^)] 2 r ^t^ 2 [ 7 | II (^)^- 1 , [t<l,R M ^oo], (4.25) 



16 ^ r 2 )l kR Ab [ k 2 r 2 [ R M 
as well as in the trigonometric region 



) 2 7^d + l) 2 , [*> 1,1 = 0, Rap ^oo], (4.26) 

^(t)--^[ln(^)] 2 ^-^, [t>l,l = 0,R Ab ^oo], (4.27) 
16 r 2 kR Ab 



PrAB(t)^ - T - A — , [t>l,l>0,R AB ^oo]. (4.28) 



3A 2 . 1 
2k 2 R 2 AB t* kR AB ' 

The supersingularity forms of P e (0, t) and P r (l, t) are obtained by first integrating 
the expressions (4.24)-(4.28) over the relevant intervals and subsequently going 
over to the limit R A b — > oo. 
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Within the exponential region, the limits t — > and Rab —> oo mutually 
strengthen the rate of vanishing. The factor (7^-)^ in the expression (4.25) van- 
ishes at B,Ab > f"2 for t — > . The functions p e AB (t) are integrable near the point 
t = even at finite values of -R^s- Moreover, its integral vanishes in the limit 
Rab — > 00. Thus 

P e As(0,£)^0, [*<l;A = a,a; J R AB ->oo]. (4.29) 

By virtue of the inequality (3.17), the condition for the absolute convergence of 
the series (3.9a) is thus fulfilled in the supersingularity limit along the exponential 
region for each of the potentials £^45 (r) of the set (4.3). 

As regards the trigonometric region, there is a competition between the 
potential and the centrifugal term and this is the point that governs convergence. In 
the absence of the latter, one extracts from Eq.(4.26) for both cases A = a, a that 

P rA/3 (l,£)^oo, [t> l;Z = O;i* A 0->oo]. (4-30) 

The convergence of the series (3.9b) is thereby frustrated for the S-wave whenever 
the physical potential decreases exponentially. Not so for the higher partial waves 
or the cases of powerlaw tails. Equations (4.27)-(4.28) imply namely that 

P rAb (l,t)^0, [t>l;l = 0;R Ab ^oo], (4.31) 

P TAB (l,t)^0, [t >l;l>0; R AB ^00}. (4.32) 

The asymptotical relationships (4.29), (4.31) and (4.32) ensure, by the the in- 
equalities (3.17) and (3.22), for the respective potential classes and partial waves, 
fulfillment of the convergence conditions (3.23), simultaneously at fixed and in- 
creasing values of the nonlinear parameter A, involved in the scattering potential 
Uab(^). One can thus write in the limit R — > 00 that 

u+(t) -> v+(t), [t < l;B = p,b;l > 0], (4.33) 
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u+(t)^v+(t), [t>l;B = b;l = 0], (4.34) 
u + (t)^>v+(t), [t>l;B = b,P;l>0]. (4.35) 

V. Asymptotical exactness 

In the last section we studied the question whether a semiclassical expansion 
that is convergent at a fixed set of dynamical parameters preserves this property 
invariably in the supersingular limit. A further point we are left with to clear is 
how the structure of the series changes upon going with the nonlinear parameter 
to infinity. A characteristic quantity of the argument for finding the answer is the 
ratio of two neighboring general terms. Recall first the definition (3.11) of q* s (t) 
along with the relevant recursion relationships (3.14) and (3.18). For simplicity, in 
the present section we are going to suppress in the formulas the matching distance 
R. Remember it is a functional of the scattering potential g 2 UA.(Rt) and is, in 
fact, invisibly present in each quantity and expression below. 

In the exponential region, the exact formula (3.14) reduces in the supersin- 
gular limit, owing to the definition (2.17), to 

«+(*)-> f dt'p^q^At'), [f<l,i*->oo]. (5.1) 
Jo 

Iteration yields then in the limit considered a familiar expression 

qtn(t)^-APe(0,t)) n , [f<l,fl->oo]. (5.2). 
n\ 

The scattering wave function reads thus by the definitions (3.9a) and (3.13) in the 
limit discussed 

vt(t) ^w+(t)expP £ (0,t), [*<l,i2->oo]. (5.3) 

For treating the trigonometric region , it is useful to introduce higher order trigono- 
metric 'coefficients', which are, as a matter of fact, no more constant. The defini- 
tion is meant to hold for both fixed or variable dynamical parameters and reads 

wt m (t) = VT (t)[C+(t)cosu T (l,t) + S+(t)smu T (l,t)}, [t> l,m = 0,l..]. (5.4) 
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Recognize the identity (5.4) reproduces at m = the definition (2.15), on account 
of which one finds that Cq = C + and Sq = S + . Combination of this identity 
with the recursion formula (3.18) furnishes in our limit, after separating the sine- 
and cosine- contributions, the following system of equations 

C+(t) - J dt'p T (t') sin W - l)]^^'), [* > 1, - oo], (5.5) 

S&(*)- J dt'p T (t')cos[kR(t' -l^w+^it'), [t>l,R^oo]. (5.6) 

Insertion of the definition (5.4) into the right hand sides of the last two formulas 
yields by analysis, owing to the infinitely rapid oscillations in the integrands, a 
pair of coupled systems of integral equations such as 

C%{t)^±£tf Pr V)S+- 1 {n [t> l,i2-oo], (5.7) 

S+(t)^-± J* tfp T {t)C+ ^(f), [t>l,R^oo]. (5.8) 

At the end of the iteration, at m = 0, one encounters the constant coefficients 
Cq , Sq . The system of integral equations becomes thereby explicitly solvable. 
The solution reads, in terms of the notation 

T m (t) = ^P T (l,t)] m , [t>l,m = 0,l,2,..], (5.9) 



as follows: 

CUt) - +T 4M (t)C+ 



S+ (f) - +T 4 ^(t)S+, (5.10) 

Sl +1 (t)^-T 4 , +l (t)C+, (5.11) 

S^ +2 (t)^-T 4tl+2 (t)S+, (5.12) 

St+3(t) - +T4„+3(t)C+, (5.13) 



where fj, = 0, 1, 2, ... Incorporation of the formulas (5.10)-(5.13) into the definition 
(3.9b) furnishes 

oo 

"r (*) "> £{PV(*) " 2W 2 (t)]«;J ) (0+ , , 

M =0 ( 5 - 14 ) 

[T^ +1 (t)-T^ +3 (t)]w; (t)}, [t > 1 }J R- oo], 
where we used the ad hoc yet not inconsistent notation 

W - Q (t) =r] T (t){S£cos[kR(t-l)}-C£sm[kR(t-l)}}, [t > 1], (5.15) 

Observe there is only a single quadrature involved in the asymptotical for- 
mula (5.14), namely the one implicitly contained in the definition (5.9). Knowledge 
of the higher order t-dependent coefficients thus rests upon the knowledge of the 
R — > oo form of the zero order constants Cq and Sq. These constants are fixed 
by the claim for the smoothness of the exact but supersingularity solution at t = 1 
as follows. 

At the matching point itself, approaching it from either of the regions e or 
r, one obtains the respective amplitude functions, see Eq.(2.16), along with the 
relevant derivatives as follows: 

[77 7 (*)]t=i->2?(fci*)3, [ 7 = e,r; R^oo], (5.16) 

[^^] t=1 ^±2HkR^[2X^-g 2 R 3 U'(R)l [ 7 = («); R^oo}. (5.17) 

The second term in the brackets vanishes in this limit by the restrictions imposed 
upon the potentials in virtue of the asymptotical relationships (2.1). 

As to the exponential region, the supersingularity formula (5.3) can be iden- 
tically recast by means of the definition (2.14) as 

v?(t) -> r) e (t)exp[uj e (l,t) +P £ (0,t)], [t<l,R^oo]. (5.18) 
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Approaching the matching point from this region , one concludes by Eqs. (5.16)- 
(5.18) that 

[v e + (*)] t =i^2*(fci2)*expP e (0,l), [R^oo], (5.19) 

[^^]t=i - 2*(fci*)* expP e (0, 1){4\1 + 2"* + [p«(f)]t=i}, [P - oo]. (5.20) 
Recall now the expressions (4.24)-(4.25) and extract from them 

[p e (*)]t=i = O{j-^}-0, [P = 6;P^oo], (5.21) 
[p e (*)]t=i = 0{-}-oo, [B = #i2-oo]. (5.22) 

P2 

Therefore, we have to restrict further discussion to scattering by potentials of the 
class UAb{f)i [A = a, a] of Eqs. (4.3). As regards the trigonometric region, the 
definition (2.15) yields 

[w+ ] t=1 - 2i(kR) 1 *C+, [B = b,R^ oo], (5.23) 

[^^}t=i^(kR^[-2^X 2 T C+ + S+], [B = b,R^oo]. (5.24) 

Smooth matching is realized by equating Eqs. (5. 19) and (5.23) as well as (5.20) 
and (5.24), respectively. In doing so, one obtains 

C+ -> expP £ (0,l), S+ ^expP e (0, 1)[2^A 2 + 2"t], [5 = b, R -> oo] . (5.25) 

As to the region e, Eqs. (4.29), (4.33) and (5.3) , as well as concerning the region r, 
Eqs. (4.31), (4.32), (5.9) and (5.14), combine to the statement that , for potentials 
UAB{f) with powerlaw tails, [B=b], the pair of zero order terms in the semiclassical 
expansions of the respective regions themselves recover in the supersingular limit 
the exact QM solution. In particular, 

u + (t) -> w+ {t), [t< 1;B = P,b; P^oo], (5.26) 

u + (t) -> w+ (t), [t > 1;B = b;l > 0];P^ oo]. (5.27) 
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In cases of exponential potential tails, the S-partial waves should be excluded from 
the present approach. 

VI. Discussion 

Usefulness of a modified semiclassical approach in treating singular scatter- 
ing has been checked in three steps. First, the conditions were reconsidered at 
which a smooth WKB method, proposed recently, pro duces convergent expansion 
of the wave function at fixed set of the potential parameters. An inherent new 
point is that, instead of a single one, a pair of integral equations should be set up, 
one for each of the exponential and the trigonometric regions. A regular solution 
working in the exponential region selects by virtue of the smoothness postulate 
a particular one out of the solutions prepared for the trigonometric region. An- 
other new item is the extension of the argument over potentials involving variable 
nonlinear parameters through the variation of which one may increase the core 
singularity to asymptotically high levels. Analyzed are four classes of interactions 
each of which is a product of a core factor implying exponential or powerlaw singu- 
larity and of a tail factor that decays either exponentially or powerlike. Indepen- 
dently of the stage of the singularity, the powerlike decaying potentials invariably 
develop absolute convergent expansions, along both regions, also in the supersin- 
gularity limit. In scattering by exponentially decaying potentials, the criteria of 
convergence seemingly fail to work within our argument of treating supersingu- 
larity limit. A further new feature of the present approach is a discussion of the 
quality of convergence. The conditions are found which the shape of the poten- 
tial has to show up so that the infinite series shrink, in the supersingularity limit, 
to a single term. It is, perhaps, worth mentioning that, by varying the length 
of the cut-off series in the region beyond the matching point, one can obtain, at 
the expense of a single quadrature, a solution that becomes correct to any pre- 
scribed order in the reciprocal supersingularity parameter at its asymptotically 
large values. 
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Recall the Born series furnishes the physical scattering wave function for 
nonsingular potentials at fixed linear and nonlinear dynamical parameters, exclu- 
sively. The proposed smooth WKB approach should work for singular potentials 
at fixed and asymptotical values of linear and nonlinear parameters involved in 
the Schrodinger equation. 

Acknowledgement Many thanks are due to Dr. G.Bencze for useful critical 
remarks. The author is grateful to Dr. G.Kluge and Dr. I. Racz for very valuable 
discussions. The work was partly supported by the Hungarian NSF under Grant 
No. OTKA 00157. 

References 

F. Calogero: Variable Phase Approach to Potential Scattering, (Academic Press, 
New York, 1967) 

F. Calogero: Phys. Rev. B 135, 693 (1964) 

G. Esposito: J.Phys A 31, 9493 (1998) 

N. Froeman and K.-F. Thylwe: J. Math. Phys 20, 1716 (1979) 

R.G. Newton: Scattering Theory of Waves and Particles, (Springer Verlag, 1981) 

T. Dolinszky: Physics Letters A 132, 69 (1988) 

T. Dolinszky: J.Math.Phys. 36, 1621 (1995) 

T. Dolinszky: J.Math.Phys. 38, 16 (1997). 



23 



